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We present measurements of the branching fractions, longitudinal polarization, and direct CP- 
violation asymmetries for the decays — >■ p^K*^ and B^ — >■ /o(980)-ft'*''' with a sample of 
(467 ± 5) X 10* BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e"'"e~ 
collider at the SLAC National Accelerator Laboratory. We observe _B+ — !> p°K*'^ with a signifi- 
cance of 5.3a and measure the branching fraction B{B+ p^K"*) = (4.6 ± 1.0 ± 0.4) x 10~^ 
the longitudinal polarization /i = 0.78 ± 0.12 ± 0.03, and the CP-violation asymmetry Acp = 
0.31 ± 0.13 ± 0.03. We observe B+ /o(980)iC*+ and measure the branching fraction B{B+ -> 
/o(980)iC*+) X S(/o(980) tt+tt^) = (4.2 ± 0.6 ± 0.3) x 10"® and the CP-violation asymmetry 
Acp= —0.15 ± 0.12 ± 0.03. The first uncertainty quoted is statistical and the second is systematic. 



PACS numbers: 13.25.Hw, ll.30.Er, 12.15.Hh 

The study of the branching fractions and angular dis- 
tributions of B meson decays to hadronic final states 
without a charm quark probes the dynamics of both the 
weak and strong interactions. It also plays an important 
role in understanding CP violation in the quark sector, 
constraining the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix parameters f?| and searching for evidence for 
physics beyond the standard model [2 Q . 

The charmless decays B pK* proceed through pen- 
guin loops and tree processes {B^ p'^K*^ is a pure 



penguin process) to two vector particles (VV). QCD fac- 
torization models predict a large longitudinal polariza- 
tion fraction Jl (of order (1 — Amy/m?^) ~ 0.9) for 
VV decays [4]. However, measurements of penguin- 
dominated VV decays give as low as ~ 0.5 Several 
attempts to understand the values of fi within or beyond 
the standard model have been made [Q]. 

For the B^ g^K*^ branching fraction, Beneke, 
Rohrer and Yang predict the CP-averaged branch- 
ing fraction to be ^ 10"*^' while Cheng and 
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Yang 3 quote (5.5lgJt^J) x IQ-^, both based on QCD 
factorization. The 90% confidence level (C.L.) upper 
limit 5+ p^K*^ branching fraction has been mea- 
sured to be < 6.1 X 10"'' [7:]. 

We report measurements of branching fractions, lon- 
gitudinal polarizations, and direct CP-violating asym- 
metries for the decay modes — > p^K*'^ and i?+ — > 
/o(980)ii:*+ where p^ and K*+ refer to the p°{nQ) and 
/'ir*+(892) resonances, respectively. The analysis is based 
on a data sample of (467 ± 5) x 10^ BB pairs, equivalent 
to an integrated luminosity of 426 fb"^ , collected with the 
BABMi detector at the PEP-II asymmetric-energy e+e" 
collider operated at the SLAC National Accelerator Lab- 
oratory. The e^e~ center-of-mass (cm.) energy is -^/s = 
10.58 GeV, corresponding to the T{AS) resonance mass 
(on-resonance data). In addition, 44.4 fb~^ of data col- 
lected 40 MeV below the T(45) resonance (off-resonance 
data) are used for background studies. We assume equal 
production rates of B^B~ and B'^B^ mesojis and charge- 
conjugate modes are implied throughout [Sj. The BaBAR 
detector is described in detail in Ref. Q. 

The B+ p°K*+ and B+ fo{980)K*+ candi- 
dates are reconstructed through the decays of or 
/o(980) ^ vr+TT-, K*+ K"^7r+ or K*+ K+tt^, with 
ifg — >■ TT+vr^ and tt" 77. The differential decay rate 
for B+ — > p'^K*^, after integrating over the angle be- 
tween the decay planes of the vector mesons, for which 
the acceptance is nearly uniform, is proportional to 

— - — sin^ 9k*+ sin^ 6pO + Jl cos^ 9k-'+ cos^ OpO , (1) 

where 9}^*+ (Opo) is the helicity angle of the K*^ (/0°)5 
defined as the angle between the daughter K (7r+) mo- 
mentum and the direction opposite to the B meson 
momentum in the K*~^ (p^) rest frame [l3|- The di- 
rect CP- violating asymmetry Acp is defined as Acp — 

(r- - r+)/(r- + r+), where r± = r(B± ^ /±) is the 

decay width for a given charged final state 

We apply the same selection criteria for p° and /o(980) 
candidates. The charged particles from the K*~^ and 
p^ decays are required to have a transverse momentum 
relative to the beam axis greater than 0.05GeV/c. The 
particles are identified as either charged pious or kaons 
by measurement of the energy loss in the tracking detec- 
tors, the number of photons recorded by the ring- imaging 
Cherenkov detector and the corresponding Cherenkov 
angle. These measurements are combined with informa- 
tion from the electromagnetic calorimeter and the instru- 
mented magnetic-flux return detector, where appropri- 
ate, to reject electrons, muons, and protons. 

The iiTg candidates are required to have a mass within 
0.01 GeV/c^ of the nominal mass [1], a decay vertex 
separated from the B meson decay vertex by at least 
20 times the uncertainty in the measurement of the sep- 
aration of the vertex positions, a flight distance in the 
direction transverse to the beam axis of at least 0.3 cm, 



and the cosine of the angle between the line joining the 
B and decay vertices and the momentum greater 
than 0.999. 

In the laboratory frame, the energy of each photon 
from the tt" candidate must be greater than 0.03 GeV, the 
7r° energy must exceed 0.25 GeV, and the reconstructed 
7r° invariant mass is required to be in the range 0.12 < 
m^j < 0.15 GeV/ c^. After selection, the 7r° candidate's 
mass is constrained to its nominal value 

We require the invariant mass of the K*^ and can- 
didates to satisfy 0.792 < rriK^ < 0.992 GeV/c^ and 
0.52 < mTT+TT- < 1.05 GeV/ c^, respectively. A B meson 
candidate is formed from the K*~^ and candidates, 
with the condition that the K*^ and p° candidates orig- 
inate from the interaction region and the of the B 
meson vertex flt is less than 100. We require that there 
is at least one additional charged track in the event and 
create a vertex for a second B meson from all remaining 
charged tracks and neutral clusters that are consistent 
with originating from the interaction region. 

The B meson candidates are characterized kinemat- 
ically by the energy difference AE = Eg — \/s/2 
and the beam energy-substituted mass rngg — 
[(s/2 + p-pb)V£;2-p|]'/', where iE,p) and 
{Eb,Pb) are the four- momenta of the T(4S') and B 
meson candidate in the laboratory frame, respectively, 
and the asterisk denotes the cm. frame. The event 
sample is taken from the region \AE\ < 0.10 GeV and 
5.255 < mEs < 5.289 GeV/c^. The extended toes range 
ensures the shape of the background distribution is 
properly modeled. Sideband events, outside the region 
\AE\ < 0.07 GeV and 5.270 < toes < 5.289 GeV/c^, are 
used to characterize the background and cross-check the 
Monte Carlo (MC) background simulations (Tlj. 

We suppress the background from B mesons decay- 
ing to charm by forming the invariant mass mu from 
combinations of two or three out of the four daugh- 
ter particles' four-momenta. The event is rejected if 
1.835 < Tod < 1.895 GeV/c^ and the charge and particle 
type of the tracks are consistent with a known decay from 
a D meson Finally, to reduce the background and 
to avoid the region where the reconstruction efficiency 
falls off rapidly for low momentum tracks, we require the 
cosines of the helicity angles of the K*^ and p° candi- 
dates to satisfy cos9ji*+ < 0.92 and | cos 0^0 1 < 0.95, 
respectively. 

To reject the background consisting of light-quark qq 
{q — u,d,s,c) continuum events, we require |cos6't| < 
0.85, where 9t is the angle, in the cm. frame, between 
the thrust axis of the B meson and that formed from 
the other tracks and neutral clusters in the event. Sig- 
nal events have an approximately uniform distribution in 
I cos 9t\, while qq continuum events peak at 1. 

After the application of the selection criteria, the aver- 
age number of p^K*^ candidates per event with K*'^ — >■ 
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-fCgTr"*" in signal MC simulations is 1.14 (1.03) for fully 
longitudinally (transversely) polarized decays. The can- 
didate with the smallest fitted decay vertex is cho- 
sen. Up to 2.1% (1.0%) of longitudinally (transversely) 
polarized MC signal events are misreconstructed, with 
one or more tracks originating from the other B me- 
son in the event. For p°K*+ with K*+ -J> K+ti°, the 
average number of candidates per event is 1.20 (1.08) 
and the fraction of misreconstructed candidates is 5.9% 
(2.7%) for fully longitudinally (transversely) polarized 
decays. For /o(980)ii'*+, the number of candidates per 
event and the fraction of misreconstructed events are 1.02 
(1.06) and 9.1% (13.8%) for decays with K*+ K°Tr+ 
{K*+ K+TT^). The and K*+ masses and widths 
in the MC simulation are taken from Ref. Isj and we use 
the measured /o(980) lineshape from Ref. |12| . 

A neural net discriminant is used to provide additional 
separation between signal and qq continuum. It is con- 
structed from six variables calculated in the cm. frame: 
the polar angles of the B meson momentum vector and 
the B meson thrust axis with respect to the beam axis, 
the angle between the B meson thrust axis and the thrust 
axis of the rest of the event, the ratio of the second- and 
zeroth-order momentum-weighted polynomial moments 
of the energy flow around the B meson thrust axis [ijj , 
the flavor of the other B meson as reported by a multi- 
variate tagging algorithm [ij], and the boost-corrected 
proper-time difference between the decays of the two 
B mesons divided by its variance. The discriminant is 
trained using MC for the signal, and qq continuum MC 
and off-resonance data for the background. 

We define an extended likelihood function to be used 
in an unbinned maximum likelihood (ML) fit as 



N 



c 



exp 



E". n 



(2) 



where the likelihood Ci for each event candidate i is the 
sum of njVj{xi; dj) over hypotheses j: two signal modes 
p^K*~^ and fo{980)K*~^ (including misreconstructed sig- 
nal candidates); qq continuum background; and nine BB 
backgrounds as discussed below. Vj [xi ; dj ) is the prod- 
uct of the probability density functions (PDFs) for hy- 
pothesis j evaluated for the i-th event's measured vari- 
ables Xi. The number of events for hypothesis j is de- 
noted by rij and TV is the total number of events in the 
sample. The quantities dj represent parameters to de- 
scribe the expected distributions of the measured vari- 
ables for each hypothesis j. Each discriminating vari- 
able Xi in the likelihood function is modeled with a PDF, 
where the parameters dj are extracted from MC simu- 
lation, off-resonance data, or (rriES, ^E) sideband data. 
The seven variables Xi used in the fit are mES, Ai5, the 
neural net output, m7r+7r-j ttikt!, the absolute cosine of 
the helicity angle of the ^"7/0(980) candidate and the 
cosine of the helicity of the candidate. Since most 



of the linear correlations among the fit variable distribu- 
tions are found to be about 1%, with a maximum of 11%, 
we take each Vj to be the product of the PDFs for the 
separate variables. 

The decays B+ ii:°7r+7r-)7r+ and B+ ^■ 

D°(^- i4r+7r°7r^)7r+ have large branching fractions and 
a similar topology to the decays under consideration. 
They are used as calibration channels. We apply the 
same selection criteria as described above except that 
the neural net is trained on the MC simulated data 
for the calibration channel under consideration; the A_E 
range is reduced to |Ai?| < 0.08 GeV; the to^^ and 
TOtt+tt- mass criteria are replaced with a mass range 
1.8445 < < 1.8845 GeV/c^; and no D meson veto 

is applied. We use the selected data to verify that the 
ML fit is performing correctly and that the MC is simu- 
lating the neural net, AE and rriES distributions. 

Backgrounds from BB decays involving charmed 
mesons are effectively suppressed by applying the veto 
on the D meson mass described above. The BB back- 
grounds that remain after the event selection criteria 
have been applied are identified and modeled using MC 
simulation. We categorize the BB backgrounds in the 
ML fit into nine main groups. Two groups represent de- 
cays where either a or a p°//o(980) is falsely recon- 
structed. Four groups represent nonresonant final states 
7r+7r-A'*+, p" iKTT)+, n+n- {K7r)+ , and /o(980) {Ktt)+, 
where (Ktt)^ stands for JsTg 7r+ or 7r°. The decays 
r]'Kl and B+ r]'K+ peak at high cos6'if.+ 
and are assigned their own category. We allow for de- 
cays from higher mass Kq{1430) states. All remaining 
BB background decays that are not accounted for by 
the above groups are assigned to a dominant remainder 
group. 

The invariant mass distributions in the ML fit are mod- 
eled with relativistic Breit-Wigner functions for the K*^ 
and /o(980), together with a polynomial of order up to 
four for the smoothly varying distribution of misrecon- 
structed candidates. Following Ref. 15|, a modified rel- 
ativistic Breit-Wigner function is used for the meson. 
The ^^0(1430) is modeled with the LASS parametriza- 
tion, which consists of the Kq{1A30) resonance together 
with an effective-range nonresonant component [l6| . For 
the signal, the distributions of the cosine of the helicity 
angles are described by Eq. [T] multiplied by a polyno- 
mial acceptance function that corrects for changes in effi- 
ciency as a function of helicity angle. The correction also 
accounts for the reduction in efficiency at helicity near 
0.78 introduced indirectly by the criteria used to veto D 
mesons. For backgrounds, the cosine of the helicity angle 
distribution is modeled with a polynomial. The neural 
net distributions are modeled using either an empirical 
nonparametric function [l7t or a histogram. For ttieSj an 
asymmetric Gaussian is used for the signal; the function 
x\/l — x-^ exp[— £(1 — x"^)] with x = rriEs/Eg and ^ a free 
parameter |18j is used for qq continuum and BB back- 
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grounds; and a combination of an asymmetric Gaussian 
with a polynomial is used for all other hypotheses. For 
AE, two Gaussians are used for signal and polynomials 
for all other hypotheses. 

We simultaneously fit for the branching fractions B, 
Acp, and /l (for B+ p°K*+ only). We allow the 
yields for all hypotheses to float except for [Ktt)'^ 
and /o(980) {Kii)'^ which are fixed to their predicted 
MC yields, assuming a branching fraction of 1 x 10^^. 
The predicted yields for the fixed modes are less than 
one event. The PDF parameters ^ for mES, the slope 
of the AE distribution, and the polynomial coefficients 
and normalizations describing the mass and helicity an- 
gle distributions are allowed to vary for the qq continuum 
and BB remainder groups. We validate the fitting pro- 
cedure and obtain the sizes of potential biases on the 
fit results by applying the fit to ensembles of simulated 
experiments using the extracted fitted yields from data. 
The observed fit biases in the MC samples are subtracted 
from the fitted yields measured in the data. 

The results of the ML fits are summarized in Ta- 
ble m where we assume a branching fraction of 100% 
for /o(980) ^ 7r+ TT". For decays with K*+ KItt+ 
{K*+ K+7:°), the event sample is 7444 (12867), with 
5959 ± 96 (10727 ± 122) fitted qq continuum events and 
1266 ± 81 (1451 ± 129) events in the BB background 
remainder group. The signal significance S is defined 
as S = y/2A ln£, where Aln£ is the change in log- 
likelihood from the maximum value to the value when 
the number of signal events is set to zero, corrected for 
systematic errors by convolving the likelihood function 
with a Gaussian distribution with a variance equal to 
the total systematic error defined below. The linear cor- 
relation coefficient between the p°K*+ and /o(980)Jsr*+ 
branching fractions is 0.25. 

Figs. [Hand [2] show the projections of the fits onto mss, 
AE, the masses, and the cosines of the helicity angles 
for decays with K*+ K'Itt+ and K*+ K+Tr°, re- 
spectively. The candidates in the figures are subject to a 
requirement on the probability ratio Vsig/ ('Psig + 'Pbkg) > 
0.9, where Vsig and Vhkg are the signal and total back- 
ground probabilities, respectively, computed without the 
use of the variable plotted. 

The systematic errors on the yields and branching frac- 
tions arise from the PDFs, fit biases, /o(980) parameters, 
interference, BB background yields, and efiiciencies. The 
PDF uncertainties are calculated by varying the PDF pa- 
rameters that are held fixed in the original fit by their 
errors, taking into account correlations. The uncertainty 
from the fit bias includes its statistical uncertainty from 
the simulated experiments and half of the correction it- 
self, added in quadrature. We allow for uncertainties in 
the /o(980) lineshape by performing a separate fit with 
the /o(980) mean and width as additional free param- 
eters. The effect of possible interference between the 
/o(980) and p° is estimated by adding the /o(980) and 
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FIG. 1: Projections of the multidimensional fit onto (a) niES, 
(b) AE, (c) TT~^ TT~ mass, (d) tt"*" mass, (e) |cos0^+^-|, 
and (f) cos6'j^o^+ for modes with K""^ Ks tt^. The points 
with error bars show the data; the solid line shows signal- 
plus-background; the green hatched area is the p^K*'^ signal; 
and the red dashed line is the /o(980)-/4'*"'" signal. 
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FIG. 2: 



Projections of the multidimensional fit onto (a) tties, 

+ _- „„„„ iA\ _0 



(b) AE, (c) tt"*" 7r~ mass, (d) K' n" mass, (ej |cos(7^+^- 
and (c) cos6'j^+^o for modes with K*^ — >■ K'^ tt". The same 
projection criteria and legend are used as in Fig. \T\ 



amplitudes together with a varying phase difference 
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TABLE I; Results for the measured B decays: signal yield Y (corrected for fit bias) and its statistical uncertainty, reconstruction 
efficiency (%), daughter branching fraction product IlBi{%) significance S (with statistical and systematic uncertainties 
included), branching fraction B, 90% C.L. upper limit (for modes with 5* < 6a), longitudinal polarization /l and CP- violating 
asymmetry Acp- 



Mode 




Y 


e(%) nBi{%) S{a) 


B (xlQ-S) 


UL (xlO-*^) 


h 


Acp 












5.3 


4.6 ±1.0 ±0.4 


6.0 


0.78 ±0.12 ±0.03 


0.31 ±0.13 ±0.03 


K*+ 




85 ± 24 


17.1 


23.1 


4.1 


4.6 ± 1.2 ± 0.5 


6.4 


0.74 ±0.13 ±0.03 


0.25 ±0.14 ± 0.03 


K*+ 


^ K+ 


67 ±31 


9.9 


33.3 


3.3 


4.4 ±2.0 ±0.5 


7.1 


0.94 ±0.27 ±0.03 


0.59 ±0.31 ±0.03 




/o(980)i^*+ 








9.0 


4.2 ±0.6 ±0.3 






-0.15 ±0.12 ±0.03 


K*+ 




69 ± 14 


17.9 


23.1 


6.0 


3.6 ±0.7 ±0.3 






-0.34 ±0.16 ±0.03 


K*+ 


^ K+ ^0 


91 ± 20 


11.3 


33.3 


6.8 


5.2 ±1.0 ±0.3 






0.14 ±0.12 ±0.03 



and using half the maximum change in the yield as an 
uncertainty. We test for the presence of a scalar /o (600) 
(or a) by adding it to our model, using the mass and 
width reported in Ref. [19]. The contribution of the BB 
backgrounds to the error is calculated by performing an 
ensemble of fits to the data where backgrounds are either 
removed from the fit (for those categories with a fitted 
number of events consistent with zero) , allowed to fioat 
(for the fixed backgrounds) or fixed to the expected num- 
ber of events calculated from MC. The error is calculated 
as half the difference between the default fit and the max- 
imum deviation seen in the ensemble of fits. Finally, the 
uncertainty on the longitudinal polarization affects the 
calculated yield efficiency. All these errors are additive 
in nature and affect the significance of the branching frac- 
tion results. We assume the sources of the uncertainties 
that contribute to the additive errors are uncorrelated 
when combined to form the overall branching fractions. 
The PDF parameter uncertainty contributes up to 0.4 
signal events to the systematic error and the fit bias be- 
tween 2.4 and 0.8 events, depending on the signal mode. 
We see no evidence for the /o(600) state. The /o(980) 
lineshape and interference account for up to 0.8 and 2.0 
events, respectively. The overall systematic error is dom- 
inated by the uncertainty in the BB backgrounds and, 
for p^K*^ , the systematic error on f^. The total addi- 
tive systematic error on the B'^ p^K*~^ signal yield is 
9.4 and 6.7 events for K*+ K° 7r+ and K*+ K+ 
7r°, respectively, and for B+ fo{980)K*+ it is 4.4 and 
1.3 events, respectively. 

Multiplicative uncertainties include reconstruction effi- 
ciency uncertainties from tracking (0.8% per track added 
linearly), charged particle identification (1.1% per track 
added linearly), tt*^ identification (3.0%), identifica- 
tion (1.0%), track multiplicity (1.0%), the number of 
BB pairs (1.1%), and MC signal statistics (0.2%). The 
total multiplicative branching fraction systematic error 
is 4.5% and 5.3% for decays with K*+ K'^Tr+ and 
_^ if+7r°, respectively. The multiplicative uncer- 
tainties for both sub-modes are correlated. The majority 



of the systematic uncertainties on and Acp cancel 
and the error is dominated by the uncertainty on the 
PDF parameters (0.02). The uncertainty due to the de- 
pendence of the reconstruction efficiency on the charge 
of the kaon is estimated from MC to be 0.005. The total 
systematic is calculated to be ±0.03 for all modes. 

In summary, we observe B'^ — > p^K*^ with a sig- 
nificance of 5.3cr. We measure the branching fraction 
B{B+ p°K*+) = (4.6 ± 1.0 ±0.4) x 10-^ the lon- 
gitudinal polarization = 0.78 ± 0.12 ± 0.03, and the 
CP-violating asymmetry Acp= 0.31 ± 0.13 ± 0.03. We 
observe P+ — > /o(980)i^*"'" and measure the branching 
fraction B{B+ fo{980)K*+) x S(/o(980) tt+tt") = 
(4.2 ± 0.6 ± 0.3) x 10"*^ and the CP-violating asymme- 
try Acp= -0.15 ±0.12 ±0.03. The B+ p"K*+ 
branching fraction is compatible with theoretical predic- 
tions 0,13]. 
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